Breaking out of the cycle: Including quiescence in cell cycle classification.

Single-cell transcriptomics has unveiled a vast landscape of cellular heterogeneity in which the cell cycle is a significant component. We trained a high-resolution cell cycle classifier (ccAFv2) using single cell RNA-seq (scRNA-seq) characterized human neural stem cells. The ccAFv2 classifies six cell cycle states (G1, Late G1, S, S/G2, G2/M, and M/Early G1) and a quiescent-like G0 state, and it incorporates a tunable parameter to filter out less certain classifications. The ccAFv2 classifier performed better than or equivalent to other state-of-the-art methods even while classifying more cell cycle states, including G0. We showcased the versatility of ccAFv2 by successfully applying it to classify cells, nuclei, and spatial transcriptomics data in humans and mice, using various normalization methods and gene identifiers. We provide methods to regress the cell cycle expression patterns out of single cell or nuclei data to uncover underlying biological signals. The classifier can be used either as an R package integrated with Seurat ( https://github.com/plaisier-lab/ccafv2_R ) or a PyPI package integrated with scanpy ( https://pypi.org/project/ccAFv2/ ). We proved that ccAFv2 has enhanced accuracy, flexibility, and adaptability across various experimental conditions, establishing ccAFv2 as a powerful tool for dissecting complex biological systems, unraveling cellular heterogeneity, and deciphering the molecular mechanisms by which proliferation and quiescence affect cellular processes.

Biomimetic nanodrug targets inflammation and suppresses YAP/TAZ to ameliorate atherosclerosis.

Atherosclerosis, a chronic inflammatory disease, is the primary cause of myocardial infarction and ischemic stroke. Recent studies have demonstrated that dysregulation of yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding domain (TAZ) contributes to plaque development, making YAP/TAZ potential therapeutic targets. However, systemic modulation of YAP/TAZ expression or activities risks serious off-target effects, limiting clinical applicability. To address the challenge, this study develops monocyte membrane-coated nanoparticles (MoNP) as a targeted delivery system for activated and inflamed endothelium lining the plaque surface. The MoNP system is used to deliver verteporfin (VP), aimed at inhibiting YAP/TAZ specifically within arterial regions prone to atherosclerosis. The results reveal that MoNP significantly enhance payload delivery to inflamed endothelial cells (EC) while avoiding phagocytic cells. When administered in mice, MoNP predominantly accumulate in intima of the atheroprone artery. MoNP-mediated delivery of VP substantially reduces YAP/TAZ expression, thereby suppressing inflammatory gene expression and macrophage infiltration in cultured EC and mouse arteries exposed to atherogenic stimuli. Importantly, this targeted VP nanodrug effectively decreases plaque development in mice without causing noticeable histopathological changes in major organs. Collectively, these findings demonstrate a lesion-targeted and pathway-specific biomimetic nanodrug, potentially leading to safer and more effective treatments for atherosclerosis.

Suppression of Ah Receptor (AhR) increases the aggressiveness of TNBC cells and 11-Cl-BBQ-activated AhR inhibits their growth.

Triple-negative breast cancer (TNBC) represents around 15% of the 2.26 million breast cancers diagnosed worldwide annually and has the worst outcome. Despite recent therapeutic advances, there remains a lack of targeted therapies for this breast cancer subtype. The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor with biological roles in regulating development, xenobiotic metabolism, cell cycle progression and cell death. AhR activation by select ligands can promote tumor suppression in multiple cancer types. AhR can negatively regulate the activity of different oncogenic signaling pathways and can directly upregulate tumor suppressor genes such as p27Kip1. To determine the role of AhR in TNBC, we generated AhR-deficient cancer cells and investigated the impact of AhR loss on TNBC cell growth phenotypes. We found that AhR-deficient MDA-MB-468 TNBC cells have increased proliferation and formed significantly more colonies compared to AhR expressing cells. These cells without AhR expression grew aggressively in vivo. To determine the molecular targets driving this phenotype, we performed transcriptomic profiling in AhR expressing and AhR knockout MDA-MB-468 cells and identified tyrosine receptor kinases, as well as other genes involved in proliferation, survival and clonogenicity that are repressed by AhR. In order to determine therapeutic targeting of AhR in TNBC, we investigated the anti-cancer effects of the novel AhR ligand 11-chloro-7H-benzimidazo[2,1-a]benzo[de]iso-quinolin-7-one (11-Cl-BBQ), which belongs to a class of high affinity, rapidly metabolized AhR ligands called benzimidazoisoquinolines (BBQs). 11-Cl-BBQ induced AhR-dependent cancer cell-selective growth inhibition and strongly inhibited colony formation in TNBC cells.

Induction of Aryl Hydrocarbon Receptor-Mediated Cancer Cell-Selective Apoptosis in Triple-Negative Breast Cancer Cells by a High-Affinity Benzimidazoisoquinoline.

Triple-negative breast cancer (TNBC) remains a disease with a paucity of targeted treatment opportunities. The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that is involved in a wide range of physiological processes, including the sensing of xenobiotics, immune function, development, and differentiation. Different small-molecule AhR ligands drive strikingly varied cellular and organismal responses. In certain cancers, AhR activation by select small molecules induces cell cycle arrest or apoptosis via activation of tumor-suppressive transcriptional programs. AhR is expressed in triple-negative breast cancers, presenting a tractable therapeutic opportunity. Here, we identify a novel ligand of the aryl hydrocarbon receptor that potently and selectively induces cell death in triple-negative breast cancer cells and TNBC stem cells via the AhR. Importantly, we found that this compound, Analog 523, exhibits minimal cytotoxicity against multiple normal human primary cells. Analog 523 represents a high-affinity AhR ligand with potential for future clinical translation as an anticancer agent.

Lesion-specific suppression of YAP/TAZ by biomimetic nanodrug ameliorates atherosclerosis development.

Atherosclerosis, characterized by the buildup of lipid-rich plaque on the vessel wall, is the primary cause of myocardial infarction and ischemic stroke. Recent studies have demonstrated that dysregulation of yes-associated protein 1 (YAP) and transcriptional coactivator with PDZ-binding domain (TAZ) contributes to plaque development, making YAP/TAZ potential therapeutic targets. However, systemic modulation of YAP/TAZ expression or activities risks serious off-target effects, limiting clinical applicability. To address the challenge, this study develops monocyte membrane-coated nanoparticles (MoNP) as a drug delivery vehicle targeting activated endothelium lining the plaque surface and utilizes MoNP to deliver verteporfin (VP), a potent YAP/TAZ inhibitor, for lesion-specific treatment of atherosclerosis. The results reveal that MoNP significantly enhance payload delivery to inflamed endothelial cells (EC) while avoiding phagocytic cells, and preferentially accumulate in atherosclerotic regions. MoNP-mediated delivery of VP substantially reduces YAP/TAZ expression, suppressing inflammatory gene expression and macrophage infiltration in cultured EC and mouse arteries exposed to atherogenic stimuli. Importantly, this lesion-targeted VP nanodrug effectively decreases plaque development in mice without causing noticeable histopathological changes in major organs. Collectively, these findings demonstrate a plaque-targeted and pathway-specific biomimetic nanodrug, potentially leading to safer and more effective treatments for atherosclerosis.

Systematic integration of protein-affecting mutations, gene fusions, and copy number alterations into a comprehensive somatic mutational profile.

Somatic mutations occur as random genetic changes in genes through protein-affecting mutations (PAMs), gene fusions, or copy number alterations (CNAs). Mutations of different types can have a similar phenotypic effect (i.e., allelic heterogeneity) and should be integrated into a unified gene mutation profile. We developed OncoMerge to fill this niche of integrating somatic mutations to capture allelic heterogeneity, assign a function to mutations, and overcome known obstacles in cancer genetics. Application of OncoMerge to TCGA Pan-Cancer Atlas increased detection of somatically mutated genes and improved the prediction of the somatic mutation role as either activating or loss of function. Using integrated somatic mutation matrices increased the power to infer gene regulatory networks and uncovered the enrichment of switch-like feedback motifs and delay-inducing feedforward loops. These studies demonstrate that OncoMerge efficiently integrates PAMs, fusions, and CNAs and strengthens downstream analyses linking somatic mutations to cancer phenotypes.

Single cell analysis reveals satellite cell heterogeneity for proinflammatory chemokine expression.

Background: The expression of proinflammatory signals at the site of muscle injury are essential for efficient tissue repair and their dysregulation can lead to inflammatory myopathies. Macrophages, neutrophils, and fibroadipogenic progenitor cells residing in the muscle are significant sources of proinflammatory cytokines and chemokines. However, the inducibility of the myogenic satellite cell population and their contribution to proinflammatory signaling is less understood. Methods: Mouse satellite cells were isolated and exposed to lipopolysaccharide (LPS) to mimic sterile skeletal muscle injury and changes in the expression of proinflammatory genes was examined by RT-qPCR and single cell RNA sequencing. Expression patterns were validated in skeletal muscle injured with cardiotoxin by RT-qPCR and immunofluorescence. Results: Satellite cells in culture were able to express Tnfa, Ccl2, and Il6, within 2 h of treatment with LPS. Single cell RNA-Seq revealed seven cell clusters representing the continuum from activation to differentiation. LPS treatment led to a heterogeneous pattern of induction of C-C and C-X-C chemokines (e.g., Ccl2, Ccl5, and Cxcl0) and cytokines (e.g., Tgfb1, Bmp2, Il18, and Il33) associated with innate immune cell recruitment and satellite cell proliferation. One cell cluster was enriched for expression of the antiviral interferon pathway genes under control conditions and LPS treatment. Activation of this pathway in satellite cells was also detectable at the site of cardiotoxin induced muscle injury. Conclusion: These data demonstrate that satellite cells respond to inflammatory signals and secrete chemokines and cytokines. Further, we identified a previously unrecognized subset of satellite cells that may act as sensors for muscle infection or injury using the antiviral interferon pathway.

Molecular subtypes of ALS are associated with differences in patient prognosis.

Amyotrophic Lateral Sclerosis (ALS) is a neurodegenerative disease with poorly understood clinical heterogeneity, underscored by significant differences in patient age at onset, symptom progression, therapeutic response, disease duration, and comorbidity presentation. We perform a patient stratification analysis to better understand the variability in ALS pathology, utilizing postmortem frontal and motor cortex transcriptomes derived from 208 patients. Building on the emerging role of transposable element (TE) expression in ALS, we consider locus-specific TEs as distinct molecular features during stratification. Here, we identify three unique molecular subtypes in this ALS cohort, with significant differences in patient survival. These results suggest independent disease mechanisms drive some of the clinical heterogeneity in ALS.

Revealing the biology behind MRI signatures in high grade glioma.

Magnetic resonance imaging (MRI) measurements are routinely collected during the treatment of high-grade gliomas (HGGs) to characterize tumor boundaries and guide surgical tumor resection. Using spatially matched MRI and transcriptomics we discovered HGG tumor biology captured by MRI measurements. We strategically overlaid the spatially matched omics characterizations onto a pre-existing transcriptional map of glioblastoma multiforme (GBM) to enhance the robustness of our analyses. We discovered that T1+C measurements, designed to capture vasculature and blood brain barrier (BBB) breakdown and subsequent contrast extravasation, also indirectly reveal immune cell infiltration. The disruption of the vasculature and BBB within the tumor creates a permissive infiltrative environment that enables the transmigration of anti-inflammatory macrophages into tumors. These relationships were validated through histology and enrichment of genes associated with immune cell transmigration and proliferation. Additionally, T2-weighted (T2W) and mean diffusivity (MD) measurements were associated with angiogenesis and validated using histology and enrichment of genes involved in neovascularization. Furthermore, we establish an unbiased approach for identifying additional linkages between MRI measurements and tumor biology in future studies, particularly with the integration of novel MRI techniques. Lastly, we illustrated how noninvasive MRI can be used to map HGG biology spatially across a tumor, and this provides a platform to develop diagnostics, prognostics, or treatment efficacy biomarkers to improve patient outcomes.

Investigating the Interactions of Glioma Stem Cells in the Perivascular Niche at Single-Cell Resolution using a Microfluidic Tumor Microenvironment Model.

The perivascular niche (PVN) is a glioblastoma tumor microenvironment (TME) that serves as a safe haven for glioma stem cells (GSCs), and acts as a reservoir that inevitably leads to tumor recurrence. Understanding cellular interactions in the PVN that drive GSC treatment resistance and stemness is crucial to develop lasting therapies for glioblastoma. The limitations of in vivo models and in vitro assays have led to critical knowledge gaps regarding the influence of various cell types in the PVN on GSCs behavior. This study developed an organotypic triculture microfluidic model as a means to recapitulate the PVN and study its impact on GSCs. This triculture platform, comprised of endothelial cells (ECs), astrocytes, and GSCs, is used to investigate GSC invasion, proliferation and stemness. Both ECs and astrocytes significantly increased invasiveness of GSCs. This study futher identified 15 ligand-receptor pairs using single-cell RNAseq with putative chemotactic mechanisms of GSCs, where the receptor is up-regulated in GSCs and the diffusible ligand is expressed in either astrocytes or ECs. Notably, the ligand-receptor pair SAA1-FPR1 is demonstrated to be involved in chemotactic invasion of GSCs toward PVN. The novel triculture platform presented herein can be used for therapeutic development and discovery of molecular mechanisms driving GSC biology.

Nonredundant, isoform-specific roles of HDAC1 in glioma stem cells.

Glioblastoma (GBM) is characterized by an aberrant yet druggable epigenetic landscape. One major family of epigenetic regulators, the histone deacetylases (HDACs), are considered promising therapeutic targets for GBM due to their repressive influences on transcription. Although HDACs share redundant functions and common substrates, the unique isoform-specific roles of different HDACs in GBM remain unclear. In neural stem cells, HDAC2 is the indispensable deacetylase to ensure normal brain development and survival in the absence of HDAC1. Surprisingly, we find that HDAC1 is the essential class I deacetylase in glioma stem cells, and its loss is not compensated for by HDAC2. Using cell-based and biochemical assays, transcriptomic analyses, and patient-derived xenograft models, we find that knockdown of HDAC1 alone has profound effects on the glioma stem cell phenotype in a p53-dependent manner. We demonstrate marked suppression in tumor growth upon targeting of HDAC1 and identify compensatory pathways that provide insights into combination therapies for GBM. Our study highlights the importance of HDAC1 in GBM and the need to develop isoform-specific drugs.

Neural G0: a quiescent-like state found in neuroepithelial-derived cells and glioma.

Single-cell RNA sequencing has emerged as a powerful tool for resolving cellular states associated with normal and maligned developmental processes. Here, we used scRNA-seq to examine the cell cycle states of expanding human neural stem cells (hNSCs). From these data, we constructed a cell cycle classifier that identifies traditional cell cycle phases and a putative quiescent-like state in neuroepithelial-derived cell types during mammalian neurogenesis and in gliomas. The Neural G0 markers are enriched with quiescent NSC genes and other neurodevelopmental markers found in non-dividing neural progenitors. Putative glioblastoma stem-like cells were significantly enriched in the Neural G0 cell population. Neural G0 cell populations and gene expression are significantly associated with less aggressive tumors and extended patient survival for gliomas. Genetic screens to identify modulators of Neural G0 revealed that knockout of genes associated with the Hippo/Yap and p53 pathways diminished Neural G0 in vitro, resulting in faster G1 transit, down-regulation of quiescence-associated markers, and loss of Neural G0 gene expression. Thus, Neural G0 represents a dynamic quiescent-like state found in neuroepithelial-derived cells and gliomas.

Gene Regulatory Network Analysis and Engineering Directs Development and Vascularization of Multilineage Human Liver Organoids.

Pluripotent stem cell (PSC)-derived organoids have emerged as novel multicellular models of human tissue development but display immature phenotypes, aberrant tissue fates, and a limited subset of cells. Here, we demonstrate that integrated analysis and engineering of gene regulatory networks (GRNs) in PSC-derived multilineage human liver organoids direct maturation and vascular morphogenesis in vitro. Overexpression of PROX1 and ATF5, combined with targeted CRISPR-based transcriptional activation of endogenous CYP3A4, reprograms tissue GRNs and improves native liver functions, such as FXR signaling, CYP3A4 enzymatic activity, and stromal cell reactivity. The engineered tissues possess superior liver identity when compared with other PSC-derived liver organoids and show the presence of hepatocyte, biliary, endothelial, and stellate-like cell populations in single-cell RNA-seq analysis. Finally, they show hepatic functions when studied in vivo. Collectively, our approach provides an experimental framework to direct organogenesis in vitro by systematically probing molecular pathways and transcriptional networks that promote tissue development.

N6-methyladenosine mRNA marking promotes selective translation of regulons required for human erythropoiesis.

Many of the regulatory features governing erythrocyte specification, maturation, and associated disorders remain enigmatic. To identify new regulators of erythropoiesis, we utilize a functional genomic screen for genes affecting expression of the erythroid marker CD235a/GYPA. Among validating hits are genes coding for the N6-methyladenosine (m6A) mRNA methyltransferase (MTase) complex, including, METTL14, METTL3, and WTAP. We demonstrate that m6A MTase activity promotes erythroid gene expression programs through selective translation of ~300 m6A marked mRNAs, including those coding for SETD histone methyltransferases, ribosomal components, and polyA RNA binding proteins. Remarkably, loss of m6A marks results in dramatic loss of H3K4me3 marks across key erythroid-specific KLF1 transcriptional targets (e.g., Heme biosynthesis genes). Further, each m6A MTase subunit and a subset of their mRNAs targets are required for human erythroid specification in primary bone-marrow derived progenitors. Thus, m6A mRNA marks promote the translation of a network of genes required for human erythropoiesis.

Monitoring changes in the healthy female metabolome across the menstrual cycle using GC × GC-TOFMS.

Urinary metabolomics offers a non-invasive means of obtaining information about the system-wide biological health of a patient. Untargeted metabolomics approaches using one-dimensional gas chromatography (GC) are limited due to the chemical complexity of urine, which poorly detects co-eluting low-abundance analytes. Metabolite detection and identification can be improved by applying comprehensive two-dimensional GC, allowing for the discovery of additional viable biomarkers of disease. In this work, we applied comprehensive two-dimensional GC coupled with time-of-flight mass spectrometry (GC × GC-TOFMS) to the analysis of urine samples collected daily across 28-days from 10 healthy female subjects for a personalized approach to female reproductive health monitoring. Through this analysis, we identified 935 unique volatile metabolites. Two statistical methods, a modified T-statistic and Wilcoxon Rank Sum, were applied to analyze differences in metabolome abundance on ovulation days as compared to non-ovulation days. Four metabolites (2-pentanone, 3-penten-2-one, carbon disulfide, acetone) were identified as statistically significant by the modified T-statistic but not the Rank Sum, after a false-discovery rate of 0.1 was set using a Benjamini-Hochberg correction. Subsequent analyses by boxplot indicated that the putative volatile metabolic biomarkers for fertility are expressed in increased or decreased abundance in urine on the day of ovulation. Individual analysis of metabolome expression across 28-days revealed some subject-specific features, which suggest a potential for long-term, personalized fertility monitoring using metabolomics.

Metadherin Is a Prognostic Apoptosis Modulator in Mesothelioma Induced via NF-κB-Mediated Signaling.

Therapies against malignant pleural mesothelioma (MPM) have yielded disappointing results, in part, because pathologic mechanisms remain obscure. In searching for rational molecular targets, we identified metadherin (MTDH), a multifunctional gene associated with several tumor types but previously unrecognized in MPM. Cox proportional hazards regression analysis delineated associations between higher MTDH expression and lower patient survival from three independent MPM cohorts (n = 349 patients). Through in vitro assays with overexpression and downregulation constructs in MPM cells, we characterized the role of MTDH. We confirmed in vivo the phenotype of altered MTDH expression in a murine xenograft model. Transcriptional regulators of MTDH were identified by chromatin immunoprecipitation. Overexpression of both MTDH mRNA (12-fold increased) and protein levels was observed in tumor tissues. MTDH stable overexpression significantly augmented proliferation, invasiveness, colony formation, chemoresistance, and an antiapoptosis phenotype, while its suppression showed opposite effects in MPM cells. Interestingly, NF-κB and c-Myc (in a feed-forward loop motif) contributed to modulating MTDH expression. Knockdown of MTDH expression profoundly retarded xenograft tumor growth. Thus, our findings support the notion that MTDH integrates upstream signals from certain transcription factors and mediates pathogenic interactions contributing to MPM traits. MTDH represents a new MPM-associated gene that can contribute to insights of MPM biology and, as such, suggest other treatment strategies.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

The Transcription Factor Nfatc2 Regulates ß-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets.

Human genome-wide association studies (GWAS) have shown that genetic variation at >130 gene loci is associated with type 2 diabetes (T2D). We asked if the expression of the candidate T2D-associated genes within these loci is regulated by a common locus in pancreatic islets. Using an obese F2 mouse intercross segregating for T2D, we show that the expression of ~40% of the T2D-associated genes is linked to a broad region on mouse chromosome (Chr) 2. As all but 9 of these genes are not physically located on Chr 2, linkage to Chr 2 suggests a genomic factor(s) located on Chr 2 regulates their expression in trans. The transcription factor Nfatc2 is physically located on Chr 2 and its expression demonstrates cis linkage; i.e., its expression maps to itself. When conditioned on the expression of Nfatc2, linkage for the T2D-associated genes was greatly diminished, supporting Nfatc2 as a driver of their expression. Plasma insulin also showed linkage to the same broad region on Chr 2. Overexpression of a constitutively active (ca) form of Nfatc2 induced ß-cell proliferation in mouse and human islets, and transcriptionally regulated more than half of the T2D-associated genes. Overexpression of either ca-Nfatc2 or ca-Nfatc1 in mouse islets enhanced insulin secretion, whereas only ca-Nfatc2 was able to promote ß-cell proliferation, suggesting distinct molecular pathways mediating insulin secretion vs. ß-cell proliferation are regulated by NFAT. Our results suggest that many of the T2D-associated genes are downstream transcriptional targets of NFAT, and may act coordinately in a pathway through which NFAT regulates ß-cell proliferation in both mouse and human islets.

MiR-155-regulated molecular network orchestrates cell fate in the innate and adaptive immune response to Mycobacterium tuberculosis.

The regulation of host-pathogen interactions during Mycobacterium tuberculosis (Mtb) infection remains unresolved. MicroRNAs (miRNAs) are important regulators of the immune system, and so we used a systems biology approach to construct an miRNA regulatory network activated in macrophages during Mtb infection. Our network comprises 77 putative miRNAs that are associated with temporal gene expression signatures in macrophages early after Mtb infection. In this study, we demonstrate a dual role for one of these regulators, miR-155. On the one hand, miR-155 maintains the survival of Mtb-infected macrophages, thereby providing a niche favoring bacterial replication; on the other hand, miR-155 promotes the survival and function of Mtb-specific T cells, enabling an effective adaptive immune response. MiR-155-induced cell survival is mediated through the SH2 domain-containing inositol 5-phosphatase 1 (SHIP1)/protein kinase B (Akt) pathway. Thus, dual regulation of the same cell survival pathway in innate and adaptive immune cells leads to vastly different outcomes with respect to bacterial containment.